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Abstract 

Results of a multi-color study of the variability of the magnetic cataclysmic vari- 
able BY Cam are presented. The observations were obtained at the Korean 1.8m and 
Ukrainian 2.6m, 1.2m and 38-cm telescopes in 2003-2005, 56 observational runs cover 189 
hours. The variations of the mean brightness in different colors are correlated with a slope 
dR/dV=1.29(4), where the number in brackets denotes the error estimates in the last dig- 
its. For individual runs, this slope is much smaller ranging from 0.98(3) to 1.24(3), with 
a mean value of 1.11(1). Near the maximum, the slope becomes smaller for some nights, 
indicating more "blue" spectral energy distribution, whereas the night-to-night variability 
has an "infrared" character. For the simultaneous UBVRI photometry, the slopes increase 
with wavelength from dU/dR=0.23(l) to dI/dR=1.18(l). Such wavelength dependence 
is opposite to that observed in non-magnetic cataclysmic variables, in an agreement to 
the model of cyclotron emission. The principal component analysis shows two (with a 
third at the limit of detection) components of variablitity with different spectral energy 
distribution, which possibly correspond to different regions of emission. The scalegram 
analysis shows a highest peak corresponding to the 200-min spin variability, its quarter 
and to the 30-min and 8-min QPOs. The amplitudes of all these components are depen- 
dent on wavelength and luminosity state. The light curves were fitted by a statistically 
optimal trigonometrical polynomial (up to 4-th order) to take into account a 4-hump 
structure. The dependences of these parameters on the phase of the beat period and 
on mean brightness are discussed. The amplitude of spin variations increases with an 
increasing wavelength and with decreasing brightness. The linear ephemeris based on 46 
mean minima for 2003-2005 is HJD 2453213.010(3) + 0.137123(3)^. The extensive tables 
of the original observations and of results of analysis are published in an electronic form. 
The nearby star GSC 4081-1562 was found to be an eclipsing red variable. 
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1 Introduction 



Cataclysmic variables are close binaries with a red dwarf filling its Roche lobe and a white 
dwarf. Accretion structure is strongly dependent on the characteristics of the flow as well as on 
the magnetic field of the compact star (see e.g. monographs by Warner [55] and Hellier [20]). 

The object H0538+608 = BY Cam was discovered as an X-ray source by Forman et al. [18]. 
Remillard et al. [44] classified it as an AM Herculis-type system. Mason et al. [34, 35] provided 
near simultaneous polarimetry and phase-resolved spectroscopy and noticed extreme variations 
in the accretion geometry. Andronov and Fuhrmann [11] suggested that these variations are due 
to a small asynchronism between the spin and orbital rotation of the magnetic white dwarfs. 
So there are slow periodic changes of the magnetic axis of the white dwarf (and corresponding 
changes of the accretion structure) in respect to the red dwarf with a period of 14. d 52 [50]. 
These changes are sometimes called "idlings" (or, alternately, " scrollings" , "slewings"). In 
classical polars, the white dwarf is synchronized with the orbital motion, but the orientation of 
the magnetic axis may undergo cyclic changes in respect to the line of centers of the rotating 
cataclysmic binary system. Such changes were called "swingings" (Andronov [2]), contrary to 
"idlings" in the asynchronous polars caused by different orbital and spin periods of the white 
dwarf. With changing orientation of the magnetic axis in respect to the secondary (donor) star, 
also changes the location and other characteristics of accretion columns above the photosphere 
of the white dwarf. 

From the cyclotron humps in the infrared spectra, Cropper et al. [17] estimated magnetic 
field strength of the white dwarf of 41MGs, which is larger than in AM Her itself. 

This puts the object between relatively numerous groups of DQ Her - type stars (or inter- 
mediate polars) with rapidly rotating white dwarf and the AM Her - type stars (or classical 
polars) which are in synchronism [39]. Another asynchronous system was Nova Cyg 1975 = 
V1500 Cyg [49], thus the desynchronization may be caused by a Nova outburst. Two more 
recent asynchronous polars are known (V1432 Aql and CI Ind). Sometimes these stars are 
called as the BY Cam - type stars (e.g. Mason [33], Honeycutt and Kafka [23]). 

According to the last "live" electronic version of the catalogue (dated 01.09.2006) by Ritter 
and Kolb [46], the number of classical, asynchronous and intermediate polars is 77, 4 and 38 
(from a total number of cataclysmic variables of 641), respectively. So asynchronous polars are 
relatively rare objects. To study BY Cam, Silber et al. [50] have organized the "Noah Project" 
with a large flow of data obtained during more than 40 nights. Mason et al. [36] have succeeded 
to explain contradictions between different periods using the model of alternately changing 
active poles of the white dwarf. In this paper, we report results from the "Noah-2" project 
initiated to study different types of variability of this interesting object. During the "Noah-2" 
project, we have obtained 56 nightly runs of observations, more than 40 nights initially planned 
for the "Noah-1" project, and more attention is paid to the color information. The extensive 
tables of individual observations and the model parameters are presented electronically. 

2 Observations and Comparison Stars 
2.1 Observations 

The observations were made in few observatories: 

• Bohyunsan Optical Astronomy Observatory (BOAO) of the Korea Astronomy Observa- 
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tory and Space Science Institute (KASI). The observations have been obtained with a 
thinned SITe 2k CCD camera attached to the 1.8m telescope. The instrumental V and 
R systems have been used in the mode of alternatively changed filters. To determine 
instrumental magnitudes of stars, the IRAF/DAOPHOT package [37] has been used. 

• Crimean Astrophysical Observatory (CrAO) with the 2.6m reflector ZTSh with the "wide 
R" filter (partially overlapping with V) equipped with a new version of the photometer- 
polarimeter, which allows to determine both linear and circular polarization. The com- 
puter program for the data reduction was introduced by Breus et al. [16]. 

• CrAO with the 1.25 m reflector AZT-11 UBVRI photometer-polarimeter [43]. 

• CrAO with the 38cm reflector K-380. Some results of this monitoring had been published 
separately [40]. However, in the present study, we have re-analyzed these data with new 
ones to make a more complete study. 

The VR brightness of the comparison stars based on the BOAO observations is shown in 
Table 1 (see a more detailed discussion below). 

The journal of observations is presented in Table 2 (available only electronically). In an 
addition to usually presented values, we have added the phase for the mid-run according to the 
"beat period" ephemeris [40] 

T phot . be atmin.(HJD) = 2453089.2473 + 14.568 • E (1) 

where the initial epoch corresponds to the minumum of mean brightness. The light curves are 
modulated with this period, thus the corresponding phase is important to compare light curves 
from different beat cycles (see below). 

Altogether 56 nights of observations were obtained, 7 from which are with changing filters 
(VR, BOAO, 1.8m, 33 h ), 1 in V (1.8m), 4 with 3-sec photometry (R, ZTSh, 2.6m), and 1 in a 
mode of simultaneous 5-color 10-sec photometry (UBVRI, AZT-11, 1.25m, 4. /l 7). The rest 47 
runs are of monochromatic CCD (R, K-380, 0.38m) photometry. 

2.2 Comparison stars 

The photometric BV standards in the field have been published by Henden and Honeycutt [22] , 
from which some were chosen as comparison star candidates by Sumner [52] . They are marked 
in Fig. 1 by "H" and "S" with a number, respectively. Henden [21] has published another (more 
extended) sequence also in BV only, with magnitudes, which are slightly different from that of 
Henden and Honeycutt [22]. For the R photometry, no published standards have been found. 

The brightness of the comparison star CO (GSC 4094:2109) was determined using the UB- 
VRI photometer-polarimeter by linking to the standard star BD+27°2120 [38]: U=10. m 20, 
B=10. m ll, V=9. m 95,R=9. m 74 I=9. m 67 [7]. The magnitudes of the comparison star C7 (GSC 
4081:280) from our measurements are U=15. m 405, B=13. m 468, V=ll. m 762, R=10. m 390, I=9. m 359. 
The BV magnitudes are in a fairly good agreement with another estimate B= 13. m 473, V=ll. m 761 
[22]. This star was used for photoelectric photometry at the 1.25m telescope in CrAO. For the 
CCD photometry obtained at the 1.8m telescope of BOAO, we have used the stars C1-C5 for 
the "artificial" comparison star [8, 26]. At the 38cm telescope of CrAO, the star Cll was 
used as the comparison star, and C8 as the check star. To determine relative brightness of 
comparison stars used at different telescopes, we have measured stars C7-C20 on 10 pairs of 
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CCD images obtained in the BOAO. As both BOAO V and R instrumental systems are close 
to the standard ones [27], no color correction was applied. 

The adopted brightness of the comparison star C5 (GSC 4081:509) is V=13. m 373 and 
R=12. m 595. Other published values are V=13. m 408±0. m 008 [22] and V=13. m 392±0. m 032 [21] 
are in a reasonable agreement with our data. This star was set as the "main" one, as it has 
the color index closest to that of the variable star. The magnitudes of Cll (GSC 4094:467) 
determined by Pavlenko [40] are V=14. m 55, R=13. m 96. These values have been used for the 
K-380 and ZTSh data. This V magnitude coincides with the USNO value and is intermediate 
between our BOAO-based estimate V=14. m 514 and the value 14. m 572 [22] and coincides with 
the USNO value. However, the magnitude in R differs from the BOAO value more significantly: 
R=13. m 68. Because of the difference between the photometric systems used, there still may be 
some systematic shifts, which are much smaller than the ~ l m amplitude. Thus we analyse the 
runs obtained using different telescopes separately. 

All original observations (HJD, magnitude) are presented electronically in Tables 3, 4 (V, 
R data obtained in the BOAO), 5 (ZTSh R data), 6 (UBVRI brightness, color and principal 
components, see section 6) and 7 (K-380 R data). 

The r.m.s. accuracy estimates for a single BOAO measurement from the "brightness-scatter" 
dependence are 0. m 004-0. m 007 (V) and 0. m 003-0. m 013 (R) for the maximum and minimum bright- 
ness, respectively. For the K380 data, the accuracy is 0. m 01-0. m 03 [40]. For 3-sec integrations 
at the 2.6m telescope, the accuracy is 0. m 02-0. m 03. 

2.3 New variable star GSC 4081:562 

The star C4 (GSC 4081:562, a 200 o = 5 h 42 m 41 s , 5 200 o = 60°52'19") was suspected to be an 
eclipsing variable in a range 13. m 27-13. m 45 in V and ll. m 43-ll. m 56 in R (based on nightly mean 
values from the BOAO data). An amplitude is small AV = 0. m 18, AR = 0. m 13, the color index 
varies from V-R=l. m 84 at the maximum to l. m 89 at the minimum and an eclipse duration is of 
few days. Despite the formal error estimates are smaller than 0. m 003, in the final publication, 
we round them to 0. m 01. The minimum occurred at JD 2453054. A very large color index B- 
V=l. m 856 argues for a late-type red component, possibly a giant. Using the magnitudes at the 
maximum and minimum, we have determined the magnitudes of the eclipsed body: V=15. m 32 
and R=13. m 82. The resulting color index V-R=l. m 50 corresponds to a hotter star. During the 
night, some photometric waves were observed with a typical peak-to-peak amplitude ~ 0. m 02, 
which may be caused by difference in extinction for stars with so different color indexes. A 
detailed study of this object will be published separately. 

3 Active vs inactive states 

The prototype polar AM Her is characterized by transitions between high and low luminosity 
states, which Hudec and Meinunger [24] called as an " active" and " inactive" ones, respectively. 
The star was mostly in the high state during the period 1928-1976, they have found 7 "ex- 
cursions" to the low state of a duration of 130-700 d . Recently Kafka and Honeycutt [25] have 
detected 10 inactive states in AM Her during 1990-2004, which seem to be more abundant 
than in previous decades. Based on their studies of magnetic and non- magnetic variables, they 
interpret relatively rare low states as resulted by different mechanisms for the systems above 
and below the period gap. For long-period systems, they are resulted by starspots at the sec- 
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ondary at the inner Lagrangian point [31], whereas for short-period systems, by the irradiation 
feedback on the secondary star [14, 29]. 

Another mechanism which is effective in a case of channelizing of the accretion flow close to 
the inner Lagrangian point, is the "magnetic valve" [1] combined with a "swinging dipole" [2], 
i.e. changing orientation of the magnetic axis of the white dwarf at a time-scale 3-11 yr. This 
is definitely not the case of BY Cam, where the spin-orbital synchronization has been not yet 
reached. However, in this case one can suggest variations of the latitude of the magnetic axis 
(and, correspondingly, of the axis of the accretion column). Similar "precession-like" variations 
may be expected in a case of non-coincidence of the spin and orbital rotational axes of the white 
dwarf [6]. Accretion variations at a few year-scale are seen in long-term brightness variations 
[15] and seasonal variations of the outburst cycle length of dwarf novae [13] and are usually 
interpreted as a solar-type magnetic cycles [45]. Such changes may also be partialy caused 
by minor variations of the orbital separation due to a third body which cause a 'trigger-like' 
amplified response in the accretion rate [10]. All these mechanisms produce quasi-perodic 
variations of the luminosity at a few year (or decade) time-scale and may act simultaneously. 

Turning now to BY Cam, the analysis of the published AFOEV (2005) data shows a wave 
with relatively sharp minima at JD 2451800 and 2453500, with a difference of ~ 1700 d and a 
full amplitude 0. m 8. Similar conclusions may be done from the AAVSO data, which are slightly 
shifted in time. 

Honeycutt and Kafka [25] suggested a 145 d periodicity from the CCD data obtained using 
the 41-cm Robotic Telescope (RoboScope) from 1991 to 2004. Unfortunately, these sources of 
data are not corrected to 198-min variability with a large (~ l m ) amplitude, even in a simplest 
form of a mean phase curve. Thus the scatter of the light curve is highly affected by this kind 
of variability of the star. 

The complete light curve from our observations in the filters V and R is shown in Fig. 2. 
One may note that the abrupt ~ l m variability of the nightly mean brightness sometimes 
occurs at a time scale of ~ l d (from our data), whereas "slow" transitions from "bright" to 
"faint" luminosity states and vice versa last ~ 16 d and ~ 35 d , respectively. Szkody et al. 
[53] detected a system at the "low" luminosity state ~ 17. m 5. Thus we may conclude that all 
our measurements correspond to a "high" state (according to usual terminology), and thus we 
will use terms "bright" and "faint" to distinguish between the runs with a mean magnitude 
difference of ~ l m . 

4 Brightness-Color Dependence 
4.1 Nightly mean values 

The Korean CCD data have been mainly obtained in a two-color (alternate switching filters 
VR) for an analysis of the color variations. To determine "quasi-simultaneous" VR data for the 
calculation of V-R (similarly to BG CMi [27] and MU Cam=lRXS J062518.2+733433 [28]), 
we have used a local cubic fit for interpolation of the brightness in V for the moments of 
observations in another filter R (and vice versa), see Andronov and Baklanov [8] for details. 
Thus the number of " quasi-simultaneous" data is slightly less than twice the minimum number 
of points (from two filters), as it is well illustrated by Table 2 (Journal of observations). One 
may note that generally the mean value of (V — R) ^ (V) — (R), if different sets of points 
are used (i.e. for (R) we used original values of the brightness, whereas for (V — R) only the 
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interpolated ones may be used). 

From 7 nights, two correspond to a "faint luminosity" state, and five - to a "bright lumi- 
nosity" state. Fig. 3 shows a dependence for the nightly mean values. It has to be noted, that 
the error estimates for the color index (V — R) is much smaller then of the (R), because the 
variations in two filters V and R are highly correlated, thus a color index has a much smaller 
amplitude (different for different nights) then a brightness in an individual filter. The high 
correlation (p = 0.97 ± 0.11) is present, with a best fit dependence: 

(V-R) = m 939(30) - 0.219(24) • ((R) - 13.98) (2) 

Here we have used the form y(x) = y + b ■ (x — x) of a linear depencence, which minimizes the 
error estimate of the zero point [56] . The amplitude of the variations of the nightly mean values 
(which mainly characterize the luminosity changes) in R is larger then in V, with value of the 
slope of dV/dR = AR/AV = 0.88(2) (or dR/dV = 1.29(4)). Sometimes in the literature, such 
a "slope" is called "gradient" . 

For small-amplitude variability seen e.g. in the intermediate polars and nova-like variables, 
the slope dV/dR may be converted to the color excess A(V — R) = —2.5 lg(dV/dR) in respect 
to the mean color index [54]. 

However, in BY Cam, the amplitude is large, thus we prefer to use slopes instead of color 
excesses for V-R and other colors. 

4.2 Nightly run values 

The slope value 1.29(4) characterizes only the mean brightness level. For faster variability, 
we have determined the slopes for the individual nights. They are listed in Table 8, which is 
available only electronically. They vary from 0.98(3) to 1.24(3), with a value 1.200(8) for all 7 
nights. The latest value is dominated by the variations of the mean brightness. To eliminate 
this strong dependence described above, we have removed corresponding nightly mean values, 
thus obtaining a significantly smaller slope of 1.105(11), which is consistent with mean slopes 
both for combined "all bright" and "all faint" runs. 

One may also suggest a minor effect of changing slope during the night, in a sense that 
AR/AV is smaller for the bright parts of the 3 h light curve, and is larger near the minimum. 
An examination of the "V — V, R — R" diagram shows that this change of slope occurs at 
V — V ~ 0. m 10. For corresponding branches of the diagram and all data, we have obtained the 
values of slope of 0.85(4) and 1.23(2). The latter value is much more closer to the one 1.29(4) for 
the mean brightness, indicating that the main source of emission (differently seen at different 
mean values of brightness) has much lower color temperature, than the one responsible for the 
3 h brightness variability. At the maximum, the color temperature becomes even larger. 

This may be illustrated also using the 4-harmonic fits to our light curves (Fig. 4). The 
"V,R" diagram shows a relatively narrow zone. However, sometimes there are "loops", as the 
variations in V and R are not precisely simultaneous. 

The "infrared" dependence of amplitudes of the spin variations on wavelength is opposite 
to a "blue" one observed in non- magnetic cataclysmic variables, e.g. TT Ari [54], where a ~ 3 h 
variability is interpreted by rotation of the accretion disk. The present result for BY Cam is 
in a qualitative agreement to the model of cyclotron emission with a maximum in the infrared: 
with an increasing mass transfer, increases the cyclotron emission, so the system becomes more 
"red" . The spin variability is due to periodic changes of the orientation of the accretion column 
in respect to the line of sight. The flux is dependent on additional parameter - the angle between 
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the column and the line of sight, also there may be eclises of the bottom part of the column. A 
quantitative modeling of the observed dependencies needs additional polarimetric and spectral 
observations, as the number of model parameters (mass of the white dwarf, accretion rate, 
orientation and thickness of accretion columns, density distribution, strength and configuration 
of the magnetic field, orbital inclination) exceeds the number of dependencies. This may be a 
subject of another paper. 

5 Spin variability 

5.1 Periods and their origin 

As the emission lines undergo significant changes with phases of both spin and beat periods, 
they are not good indicators of the orbital motion. A linear regression to the timings of inferior 
conjunction of the red dwarf gives 

T S p(HJD) = 2451137.3952(36) + 0.13975284(14) • E (3) 

[48], known as the "201-minute" (or "orbital") period. This differs from previously published 
estimates. The spectroscopic period is distinctly different from the photometric one with an 
ephemeris for 1991-2004: 

T phot .min.(HJD) = 2448486.98(3) + 0.137120(2) • E (4) 

[23]. This difference was explained by Mason et al. [36] as a n 2uo — Q" frequency, where a; is a 
spin frequency of the white dwarf, and Q is an orbital frequency, following the model by Wynn 
and King [57] for intermediate polars. This photometric period coincides with that from the 
more recent ephemeris 

T phot , min {HJD) = 2453022.268 + 0.137120 • E (5) 

(Pavlenko 2006). Using all available timings for the main maximum, Pavlenko et al. [41] have 
published an ephemeris 

T spin (HJD) = 2453089.307 + 0.1384274 • E. (6) 

We also will use the "beat" ephemeris correponding to a synodic (spin - orbital beat) period 
Pbeat = {Pspln ~ Pml) 1 (Eq.l). These two ephemerids have been used for further study. 

5.2 4-harmonic model 

For periodic light curves, usually a multi-harmonic fit 

m{t) = d + j^{C 2j sm(2njt/P) + C 2j+1 cos(2njt/P)). (7) 

3=1 

is applied with a determination of statistically optimal value of the degree of the trigonometric 
polynomial s. Here m(t) is a signal at the moment t, P is an adopted period (we have used 
P = 0. d 137120 published by Pavlenko [40]), and Ck, k = l..(2s + 1) are parameters determined 
using the least squares method. A description of working formulae for the computer program 
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was presented by Andronov [3]. Suggesting a 4- hump model of the light curve [40], which is 
often seen in our data, we have applied the fixed value of s = 4. 

The individual light curves for different telescopes are shown in Fig. 5-7 together with the 
corresponding fits. The characteristics of the maxima and minima will be discussed in section 
7. They are presented in the electronic Table 9. 

In Fig. 5, the fits are shown for the V, R and V-R curves obtained in BOAO. They show 
significant variability of the shape, which makes senseles to put all the data simply to an 
ordinary joint phase curve. However, the curves at the "bright" state show a smaller scatter 
as compared to two curves in the "faint" state. These two curves are upper ones at the V-R 
diagram because of the "mean brightness - color" dependence mentioned above. 

6 UBVRI Photometry 

6.1 Color gradients 

The light curve obtained on January 17, 2004 at the AZT-11 telescope of CrAO is shown in the 
upper part of Fig. 8. It shows two unequal minima, onto which are superimposed possible quasi- 
periodic oscillations and flickering. At the deeper minimum, the photon counting statistics is 
much worse, which causes an increase of observational errors. 

The color indexes (middle part of the Fig. 8) show that at minima, the emission becomes 
"bluer" in a sense of decrease of the color indexes. The amplitude of orbital variations in- 
creases with wavelength. The slopes are dU/dR=0.23(l), dB/dR=0.43(l), dV/dR=0.78(l), 
dI/dR=1.18(l). This is opposite to the behavior of eclipsing cataclysmic variables with an 
accretion disk [55]. Qualitatively this situation may be explained by a "redder" spectral energy 
distribution of the accretion column than that of the white dwarf. 

6.2 Principal Component (PC) Analysis (A) 

Andronov et al. [12] have applied the PCA to determination of statistically significant number 
of sources of variability and their characteristics in 3 cataclysmic variables with different degree 
of influence of the magnetic field onto accretion. This method is an effective tool to determine 
components of variability, which have different spectral energy distributions, i.e. the wavelength 
dependencies of the amplitudes of variability. For two-color observations of the same accuracy 
and subtracted mean value, this method is known as the "orthogonal regression"; the principal 
components Uk a are Uki = {xki + Xk2)/V^, U^t. = (xki — x k2)/V2, where Xk a are observations at 
point k in channel (filter) a. So the first principal component shows common features in different 
channels (although with channel-dependent amplitude), and other PCs describe differences. 
The "color-brightness" diagrams, which are often used in astronomy (the Herzsprung-Russel 
diagram, diagrams for variable stars), are also a modification of the PCA with an aim to 
get more information on the objects. Polarimetry is also related to PCA, being applied to 
"ordinary" and "extraordinary" waves in the emission. 

For the multi-color photometry, PCA is a supplementary statistically optimal tool to extract 
components of variability, which may be overseen if analyzing the colors only. Formerly, a 
method was proposed [54], which allowed to extract only the first component. 

Applying the formalism of PCA [12] for the present data, we have found that the difference 
in the amplitude of the first principal component (Ul) is much larger than of the others. It 
explains 84% of overall variance in all colors. The second component U2 is responsible for 



8 



10%, so the total contribution of the rest three (U3, U4, U5) is only 5.8%. Despite a scatter 
of the third component, it has obvious variations at the orbital time scale. The amplitude 
"signal/noise" ratio obtained using the "running parabola" scalegram analysis [4] is 26, 14, 8, 4 
and 4 for U1-U5, thus we conclude that first 3 principal components are statistically significant. 
Relatively large values for U4 and U5 are due to larger errors at a second more deep minimum. 

The dependence of the matrix Vij on the number of channel (filter) i and on the number of the 
PC j is shown in Fig. 9 and is listed in the table below. The first (largest in the amplitude) PC 
shows both 3-hour and faster variability, as one may see in Fig. 8. Its contribution monotonically 
decreases with the number of channel (and corresponding wavelength), as one may see from 
Fig.8. The ratios Vu/V^i (which are equal to ratios of the amplitudes of variability with a shape 
Uki are equal to 0.29, 0.47, 0.83, 1 and 1.20. The difference of these values from the ratios 
dU/dR, dl/dR is present, because 3 components are statistically significant (for an exact 
linear dependence of the brightness variations in different colors, i.e. the presence of only one 
statistically significant component, the values in a pair for a given color must coincide). 

The first PC Uki mostly resembles variations in I and R and much worse the variations in 
U, what is expected, looking on the increase of Vn with wavelength. Such a behaviour may 
be interpreted for this system, assuming that in BY Cam there are sources of emission with 
different spectral energy distribution: the dominating source of variability (first PCA) is due 
to cyclotron emission of the accretion column, so it has an amplitude, which strongly increases 
with wavelength. This resembles a classical polar AM Her with an accretion column ("infrared- 
dominated" variability), but quite opposite to the nova-like variable TT Ari ("ultraviolet- 
dominated" variability) with an accretion disk. 

For the second PC in BY Cam, the dependence of Va on i has an opposite sign of slope 
than Vn, i.e. V i2 monotonically decreases with wavelength. This means that the second (in 
the amplitude) component has a "blue" spectral energy distribution, contrary to the "infrared" 
first PC. An examination of the curve (Fig. 8) shows that the second PC Uki is dominated by 
~ 30— min variability, whereas Uki is "spin-dominated". 

6.3 Scalegram Analysis 

The " A— scalegrams" [5] are shown in Fig. 10 for the original UBVRI data, for the corresponding 
color indexes and for the PCs. All the curves have the same scale, but shifted to avoid an 
overlap. This method is some kind of the wavelet analysis, which is most effective for quasi- 
periodic oscillations, where one expects to determine a characteristic cycle length, rather than 
a constant period. Contrary to the periodogram, wavelet and auto-correlation methods, which 
are higly affected by slow (periodic or aperiodic) trends, the "A— scalegram" is insensitive to 
polynomial trends up to degree 3. Thus, in a presence of strong asinusoidal spin variability (as 
in BY Cam and many other cataclysmic variables), we prefer to use this method to estimate 
effective (mean) amplitudes and periods (or cycles) of different components of variability. 

The height of the peak is proportional to the square of the "effective amplitude". The most 
pronounced peaks are seen for the scalegram for the first PC (Ul) corresponding to "effective 
periods" of 0. d 122, 0. d 37, 0. d 21 (in order of decreasing height at the scalegram). The positions 
of the peaks are shown by vertical lines for comparison with that for other scalegrams. The 
highest peak corresponds to the most prominent "3- hour" variability. It is slightly less than a 
photometric period, what is expected for a non-sinusoidal shape of the curve [5]. The second 
period is close P p hot/^, confirming a 4-peak structure of the light curve, which was pointed out 
by Pavlenko [40]. 
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The third peak corresponds to ~ 30— min QPOs. Previously, QPOs in BY Cam were 
reported to occur at some phases with different estimates of their amplitude and time scale: 5 
min [44], 10 and 15 min [50], 5-10 and 30 min [42, 51]. The fourth peak occurs at 0. d 0056=8 min 
and agrees with smaller periods mentioned above. However, its amplitude is smaller than for 
other peaks. Despite the peaks seem too small as compared to the main peak, the corresponding 
effective semi-amplitudes are large enough (0. m 58, 0. m 33, 0. m 24 and 0. m ll for four peaks for Ul, 
respectively). 

For U2, the peaks are of much smaller amplitude, but nearly at the same places, as for Ul. 
The highest peak of A (At) corresponds to 30-min component, similar to the scalegram for the 
filter U. This also is in an agreement with the conclusion of the "infrared" Ul and "blue" U2. 
The third U3 has a small amplitude, however, its variations are mainly concentrated to the 
"3- hour" peak. 

In Fig. 10, the scalegram for 3-sec photometry (all 4 nights) obtained at 2.6m telescope 
ZTSh is also shown. Besides the dominating peak, it also has smaller peaks close to the same 
positions. The " A— scalegrams" for other runs show night-to night variability of the shape and 
corresponding amplitudes and "periods". This causes an overlap of the peaks and decrease of 
the ratio "peak-to-continuum". Thus an absence of well pronounced secondary peaks agrees 
with the former conclusion of significant variability of the characteristics of these brightness 
oscillations either during the run, or from run to run. 

The height of the primary peak, which corresponds to the 3 h variability, characterizes the 
amplitude of variations. Using the 4-harmonic fits, this value was determined for all runs, which 
sufficiently cover the photometric phase, and will be discussed in the last section in more detail. 
We have splitted the "all" data sets for the BOAO and K-380 telescopes into the "bright" and 
"faint" parts. 

For the BOAO data, the height of the peak and a corresponding effective amplitude is larger 
for the filter R than for the filter V for "all", "bright" and "faint" data sets, in an agreement 
with the direct determination of the slope dV/dR< 1. 

The effective amplitude is by a factor of ~ 1.3 larger at the "faint" state than at the "bright" 
one for both BOAO and K-380 data. The position of the peak in the "faint" state switches 
to the half of the period, indicating for the redistribution of the amplitudes of the harmonic 
components of the light curve. 

7 Characteristics of Individual Light Curves 

As one of the main results of the first "Noah project", Silber et al. [50] pointed out, that the 
phase of the sine wave fits for the 3 h variability show an abrupt jump with a period of 7. d 26, 
which was interpreted as a half of the complete beat period [36]. So there are two switchings of 
the accreting column from one pole to another. If the 7. d 26 - period dominates, this means that 
two accreting poles at the white dwarf are separated by ~ 180° in longitude (with a similar 
latitude) and there will be two waves during one beat period. Diversity of the location or 
structure will lead to an unequality of these waves. 

Pavlenko [40], based on 31 runs in the inactive state, has suggested a more complicated 
structure with 3 possible jumps at the phases fbeat = 0.25, 0.55 and 0.78 with a possible interval 
of the "erratic" phase variations between the beat phases 0.0 and 0.25. The mean slopes of 
the O-C diagrams between the "main" switches at tpbeat = 0.25, 0.78 correspond to the "true" 
spin period P sp i n = 0. d 138428. However, an additional jump at (fbeat = 0.55 was interpreted as a 
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switch between 3 (or even 4) accretion zones. In the beat phase intervals 0.25-0.55 and 0.55-0.78 
for BY Cam, Pavlenko [40] had suggested a larger slope "phase/beat phase" dip/dipb ea t) (and 
thus an apparent photometric period, as P = Pq(1 + (Po/Pb ea t) dtp /dip beat), where P is a trial 
value of photometric period, Pq— its adopted initial value, Pbeat— beat period). 

It should be noted that in some nearly-synchronous polars the accreting region moves rather 
smoothly in longitude even while accreting onto a single pole, as the treading region drifts 
due the asynchronism. At first, such a model was discussed by Geckeler and Staubert [19] 
for another asynchronous polar V1432 Aql (see Andronov et al. [9] for results of the recent 
observational campaign on that object). Periodic variations of phases of the spin variability 
with the orbital phase were reported by Kim et al. [27] for another magnetic cataclysmic 
variable 1RXS J062518. 2+733433 = MU Cam. The system is an intermediate polar, and such 
a dependence may argue for a diskless model of the system [32]. Despite the spin period in BY 
Cam is close to the orbital one, contrary to the intermediate polars, the asynchronous polars 
may be also called "nearly synchronous" intermediate polars. So asynchronous polars (or BY 
Cam-stars [33,23]) are "intermediate" between the classical and intermediate polars. The ratio 
of the spin period to the orbital one for the magnetic white dwarfs in cataclysmic binaries leads 
to a phenomenological classification [39], but it has a strong physical explanation [20,32,55]: for 
small (for a given accretion rate and orbital separation) magnetic moments, the thread radius 
is small, and there is a strong accretion disc without a strong evidence for accretion columns 
("non- magnetic" systems, e.g. TT Ari). With increasing magnetic moment, increases the inner 
radius of the accretion disc, and accretion columns become stronger ("disc-fed" intermediate 
polars). Accretion disc may be disrupted in "disc- less" intermediate polars, but the stable 
spin period is shorter than the orbital one. Asynchronous polars evolve to a "phase- locked" 
synchronous rotation, which is typical for classical polars of AM Her - type. However, because 
of influence of other system parameters, sometimes systems with larger magnetic field of the 
white dwarf may stay at the "less-magnetic" branch of the classificatin described above. 

The main characteristics of the fits for BY Cam have been plotted as a function of the 
beat phase (Eq. 1) in Fig. 11. The statistically significant degree s of the polynomial fit was 
determined for each run separately, and ranges from s = 1 (sinusoidal variation) to s = 4. They 
and the light curve at the 10 equidistant phases are presented electronically in Table 9. Each 
point represent one night of the observations. As the majority of data were obtained in the 
photometric system R, we have used this sample for better statistics. Incomplete curves, where 
the data cover less than 0. d l have been excluded from the analysis, to avoid large errors of the 
fit in gaps. The resulting number of nights is 46. 

To reduce the effect of inhomogeneous distribution of phases of observations, we have used 
the "orbit-averaged" mean m mean (i.e. the coefficient C\ in the fit) instead of the "sample 
mean" . 

As one may see from the dependence of this "mean" value, the "bright" and "faint" states 
are well separated from each other. There are only 3 nights (from 45) inbetween these two 
groups. The transitions are much less frequent (i.e. 6%) than the bright (n=ll, 24%) and 
"faint" (n=32; 70%) states. The range of the mean values in the bright state (~ 0. m 3) is much 
smaller than in the faint state (~ l. m l). This is also valid separately for the maxima ("max" 
in Fig. 11) and minima ("min"), arguing for a better stability of the light curve in the bright 
state. 

Not unexpectedly, the maxima and minima in the bright state are brighter than these 
characteristics for the runs in the faint state. In the faint state, the "mean", "max" and "min" 
values show minima at the beat phases <Pbeat = and 0.5. The initial epoch was chosen by 
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Pavlenko [40] to correspond to the minimum of the sample mean brightness, so our result for 
the "orbital mean" just shows that the difference between the "sample" and "orbital" mean 
values is much smaller than the amplitude of the curve. 

The "max" and "mean" values m min and m max are correlated (p = 0.91(6)) with a sample 
mean values of 13. m 97(3) and 15. m 04(3). If (as for truly monoperiodic stars) there were no 
phase shifts, these values could be named as the " maximum" and " minimum" value of the mean 
phase curve. The slope dm max /dm m i n =0.74(5) is significantly smaller than unity, indicating 
an increase of the amplitude {m min — m max ) on m min . More usual form is the "amplitude 
) vs. m mean " . The correlation coefficient p = 0.40(14) is much smaller, but the 
correlation is still significant at the "3<r" level. The corresponding mean values are 14. m 53(7) 
(13 m 59-15 m 46) and l. m 07(3) (0 m 74-l m 46) for m mean and (m min - m max ), respectively. The 
slope is positive d(m min - m max ) / dm mean = +0.17(6). 

To remove an obvious correlations between the brightness at some phases owed to variations 
of the mean brightness and to check the stability of the light curve, we have introduced two 
related parameters Tn max n^mean and in-min m mean ■ 

For this further analysis, we have classified "intermediate" points as "bright", thus deter- 
mining a dividing value of the mean brightness of R=14. m 21. 

The parameter m max — m mean (from -0. m 6 to -0. m 3) for the bright state is much closer to 
zero than that for the faint state (from -l. m to -0. m 4). This indicates a more extended hump 
for the bright state and more stable shape of the light curve. For the (m min — m mean ) values, 
the range (0 m 4-0 m 7) for the "bright" state is also smaller than for the "faint" (0 m 29-0 m 88), 
but there is no significant systematic shift. 

Both characteristics [rn max — m mean ) and (m min — m mean ) have a large scatter. In the beat 
phase interval 0.03-0.15, there are 5 outstanding points (nights) with an amplitude, which is 
nearly twice larger than that at the nearby beat phases. The mean values for 46 runs are 
{ijn max - m mean )) = -0 m 57 (from -0 m 83 to -0 m 31) and {{m min - m mean )) = +0. m 52 (from 
0. m 29 to 0. m 88). Unexpectedly, there is no correlation between these two characteristics. If one 
will define fh mean = (m max + m min )/2, then the values (m max - fn mean ) = -{m min - m mea „), 
and there will be an exact negative correlation. The absence of correlation argues for a high 
variability of the shape of the light curve. 

The dependence of the photometric phases ip p h on the beat phases (fibeat is shown at the 
bottom part of Fig. 11. In the previous studies, only one phase of the maximum was taken 
into account (one-harmonic fit by Silber et al. [50] or an "asymptotic parabola" fit by Pavlenko 
[40]). Assuming a complicated structure of the light curve and applying a 4-harmonic fit, we 
have used the phases of both maximum <p m i n and minimum ip max and the phase difference 
between them. The last value is called "asymmetry" and is an important parameter listed in 
the General Catalogue of Variable Stars [30, 47]. 

The corresponding dependencies have following differences. The parameter <p max has jumps 
at (fbeat = 0.78 and 0.23, as was found by Pavlenko [40]. However, all the parameters show out- 
standing points at corresponding diagrams, which are characterized by relatively small formal 
error estimates, and thus may be classified as statistically significant. The original brightness 
variations, which cause such pecularities, may be produced by cloud-type inhomogeneties of 
the accretion stream. Such "flare-like" events with an amplitude of hundreds millimagnitudes 
are distinctly seen at our curves. So a "filtration" of the "bad" curves producing outstanding 
points may be done, and more smooth curves could be seen the diagrams (Fig. 11). The 
criteria for such a "filtration" may be elaborated in future. In this work, we present results of 
modeling of all curves in extensive electronic tables as "observational facts", not removing any 
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outstanding data. 

In other words, present results show that the situation is much more complicated, not 
because of bad observations, but because of the complexity of the system. 

At present, we may conclude that outstanding points really exist, but, if removing them, 
one may suggest good sequences for the remaining points. One can naively expect, that the 
light curve is dependent purely on the beat phase (thus on the periodically varying orientation 
of the magnetic field and the structure of the accretion flow). The scatter at the diagrams 
(Fig. 11) argues for additional valuable processes, particularly, mean magnitude changes, and 
large-amplitude flares. 

For our data, we suggest that the phases of minimum show the main jump at slightly larger 
phase upbeat = 0.92 for the "faint" data and possibly a more smooth variation for the "bright" 
data. There are obvious outstanding points, which correspond to the "global" lowest minimum, 
whereas at the phase light curve there may occur two local minima inside a "main" minimum. 
This does not challenge a general conclusion of switching accretion from one pole to another 
while " idling" of the white dwarf in respect to the secondary. 

The mean value of ^ min is 0.058(24), so the deviation from zero is formally not statistically 
significant. However, the initial epoch in Eq. (5) may be corrected to T p h ot min (H J D) = 
2453022.276(3) (assuming the period of Pavlenko [40]). Our 46 minima timings are best fitted 
to the ephemeris 

T phot . mm {HJD) = 2453213.010(3) + 0.137123(3) • E. (8) 

So the period value differs from that of Pavlenko [40] by la, i.e. not significantly. The maximum 
phase shift due to the difference between the periods for our 446 d of observations is only 0.01(1). 
This is comparable with the mean accuracy of the phase determination for a single curve. Thus 
a recomputation of all phases is not necessary. 

The asymmetry ip max — (Pmin has a drastic range of variations from 0.18 to 0.80, which is 
much larger than a mean accuracy estimate of 0.012. There is a slight evidence (at the 3a 
limit) for a double-hump structure with minima close to the phases of jumps. 

8 Results 

• The main differences between the spin light curves in the "bright" and "faint" brightness 
states are: a) the mean color index V-R is larger for the "bright" state (this argues that 
the emission region associated with accretion has an "infrared" energy distribution); b) 
the "bright" light curve has a smaller amplitude; c) the "bright" light curve is more 
stable; d) the "bright" light curve has a smaller value of \m max — m mean \ indicating larger 
duration of the maximum (which is usually split into two humps). 

• the slopes dV/dR are more close to unity for the variability during a night, than for 
the mean brightness, indicating that the variability due to changing orientation is more 
"blue" than the emission from the mean accretion component and is more "red" than 
the total emission. This may be tentatively interpreted as a sequence with decreasing 
temperature "white dwarf +red dwarf, "accretion structure variable with orientation", 
and "accretion structure, which is not variable with orientation". However, at present, 
this conclusion is far from the theoretical modeling. 
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• Contrary to the intermediate polars BG CMi and MU Cam, the variations in different 
colors show an increase of the amplitude with wavelength. 

• The second principal component of variability is characterized by a decrease of the am- 
plitude with wavelength. The amplitudes for the 30-min, P/4 and P variations are 
comparable, indicating that 30-min QPOs are more "blue" than the spin variability. 

• The characteristics of 46 pairs of spin maxima and minima for the season in 2003-2005 
have been determined. 

• The optimal mathematical model for the four-hump light variations corresponds to a 
fourth-order trigonometrical polynomial. The dependences of these parameters on the 
phase of the beat period and on mean brightness are discussed. In an addition, the 
30-min and 8-min quasi-periodic oscillations are present with a variable amplitude and 
effective cycle length. 

• The ephemeris for the orbital minima has been corrected. 

The determined parameters may be used for the comparison to theoretic models of accretion 
in magnetic cataclysmic variables. The further multi-color monitoring is needed for studies of 
variability of the changes of the accretion structure and rotation of the white dwarf. 
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Figure 1: Finding chart for BY Cam. The size of the field is 5.6', North is up and East is left. 
For a cross-identification, the stars from Henden and Honeycutt [22] (marked with " H" 

and Sumner [52] (marked with " S" ) are shown. 
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Figure 2: Dependence of nightly mean brightness in the filters V and R on time. R observations 
are taken for all 4 telescopes, whereas V are shown only for the 1.8m BOAO telescope. The 
error estimates of the correponding color index V-R is much smaller than of the original colors 
because of much smaller amplitude of its variability. The transitions between the "bright" and 
"faint" states usually take weeks, but sometimes even one day. 
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Figure 3: Dependence of the nightly mean color V-R on the brightness in R for the BOAO 
observations and the best fit (2). 
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Figure 4: The "V-R" diagram for the 4-harmonic fits to the BOAO data. The phase distance 
between two subsequent points is 0.01P=2 min. The line R — R = V — V with an unitary slope 
is shown for reference of individual slopes within 7 nights. 
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Figure 5: The individual light (V,R) an d color index (V-R) curves (points) for the BOAO data. 
The original data are plotted once and are not repeated once per unit phase, but the fit is 
repeated for clarity. The four-harmonic fits for individual runs are shown by solid lines. Strong 
night-to-night variability of the shape is present, with an easy visual separation of two curves 
observed at the "faint" state and five at the "bright state". 
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Figure 6: The individual light (R) curves (points) for the data obtained at 2.6m Shain telescope. 
As in Fig. 5, The original data are plotted once and are not repeated once per unit phase, but 
the fit is repeated for clarity. The four-harmonic fits for individual runs are shown by solid 
lines. 



Table 1: Brightness of BY Cam and the comparison stars from the BOAO observations (VR). 
The values of o correspond to the r.m.s. scatter of a single observation in respect to the 
"artificial comparison star". 
Star V ay R or 
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0.001 
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0.001 
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0.007 
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0.007 
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0.017 
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0.005 
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Figure 7: The individual light (R) curves (points) for the data obtained at the telescope K-380. 
As in Fig. 5, The original data are plotted once and are not repeated once per unit phase, but 
the fit is repeated for clarity. The four-harmonic fits for individual runs are shown by solid 
lines. 
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Figure 8: The individual light (UBVRI), color index (U-B, B-V, V-R, R-I, B-R) and principal 
component (U1,...,U5) curves (points) for the data obtained at the telescope AZT-11. The 
direction and scale are the same for all graphs. The interval between minor tics is 0. m 2. 
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Figure 9: Dependence of the matrix of eigenvectors V^- on the number of channel i and on the 
number of the principal component j (see Andronov et al. [12] for a detailed description). 
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Figure 10: "A"— scalegrams [5] for the individual light (UBVRI), color index (U-B, B-V, V-R, 
R-I, B-R) and principal component (U1,...,U5) curves for the data obtained at the telescope 
AZT-11 and (bottom) for all 4 nights obtained at ZTSh. The vertical lines mark positions of 
3 most prominent peaks at the scalegram for the first principal component Ul. 
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Figure 11: Dependence of the main parameters of the 4-harmonic fits to individual runs on 
the phase of the beat (synodic) period. Each point represent one night of observations. For 
the phase range 0-1, the data are marked by asterics and crosses for the "bright" and "faint" 
states, respectively. The error bars are shown, but at the left part of the graphs, they are often 
smaller than the symbols, so at the right part (phases from 1 to 2) only error bars are shown. 
The slope of the inclined lines correspond to the spin period of the white dwarf, the jumps are 
interpreted as switchings of accretion from one pole to another. 
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Table 2: Journal of observations of BY Cam: JD- integer part of the starting Julian date 
(used for a legend of the run, e.g. 52634V); HJD st(lrt - begin of observations (HJD-2400000); 
telescope (Z=2.6m, K=1.8m, A=1.25m, P=0.38m); filter; exp - exposure in seconds; N- number 
of observations; duration of the run in days. In the electronic version of the table, there 
are additional columns: HJD fi n i S h end of observations (HJD-2400000); magnitude range for 
individual data points m max , m min ; nightly mean (m) and it's accuracy estimate <r((m)); r.m.s. 
deviation of the single observation from the mean a(m); ip b - phase of the run center according 
to the "beat period" ephemeris T = 2453089.2473 + 14.568E [40]. 
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Table 2 (continued). 
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Table 3: Table of V observations (HJD-2400000, magnitude) obtained at the 1.8m telescope 
of the BOAO, Korea. The table is very long, and is planed to be published as a text file 
electronically only 

Table 4: Table of R observations (HJD-2400000, magnitude) obtained at the 1.8m telescope 
of the BOAO, Korea. The table is very long, and is planed to be published as a text file 
electronically only 

Table 5: Table of R observations (HJD-2400000, magnitude) obtained at the 2.6m telescope 
ZTSh of the CrAO, Ukraine. The table is very long, and is planed to be published as a text file 
electronically only 

Table 6: Table of brightness, color and principal components (HJD-2453022, U,B,V,R,I, U- 
B,B-V,V-R,R-I,B-R,U1,U2,U3,U4,U5) obtained at the 1.25m telescope AZT-11 of the CrAO, 
Ukraine. The table is very long, and is planed to be published as a text file electronically only 

Table 7: Table of R observations (HJD-2400000, magnitude) obtained at the 0.38m telescope 
K-380 of the CrAO, Ukraine. The table is very long, and is planed to be published as a text file 
electronically only 



Table 8: Color gradients dR/dV and dV/dR, number of pairs of interpolated VR measurements 
n, mean magnitude in V and R, correlation coefficient and the ratio to its statistical error (for 
individual nights, bright and faint state and data with a substructed corresponding nightly 
mean). The table is planed to be published as a text file electronically only. 



Table 9: The characteristics of the statistically optimal (up to 4-harmonic) fits to the individual 
runs and their error estimates: ipbeat, t mean , m mean , number of the statistically significant degree 
of the trigonometric polynomial s, unbiased estimate of the "unit weight error" s , r.m.s. 

accuracy Of the fit Sc fmaxi ^maxi T^max f^meant tfnaxi ^Pmini ^mini ^min ^mini tmini ^Pmax 

<^min, rn m i n — m max ; the light curve m(ip) for if = 0.0, 0.1,... 0.9. The table is long and wide, 
and is planed to be published as a text file electronically only. 
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